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We present an efficient method of imaging 3D nanoscale lattice behavior and strain fields in crystalline mate¬ 
rials with a new methodology - three dimensional Bragg projection ptychography (3DBPP). In this method, the 
2D sample structure information encoded in a coherent high-angle Bragg peak measured at a fixed angle is com¬ 
bined with the real-space scanning probe positions to reconstruct the 3D sample structure. This work introduces 
an entirely new means of three dimensional structural imaging of nanoscale materials and eliminates the ex¬ 
perimental complexities associated with rotating nanoscale samples. We present the framework for the method 
and demonstrate our approach with a numerical demonstration, an analytical derivation, and an experimental 
reconsh'uction of lattice distortions in a component of a nanoelectronic prototype device. 


Inversion methods provide a powerful alternative to tradi¬ 
tional objective-lens-based microscopy. Techniques that nu¬ 
merically invert coherent diffraction patterns into real space 
images have provided substantial gains in resolution and sen¬ 
sitivity in certain optical, electron, and x-ray microscopy 
experiments, especially where image-forming lenses are in¬ 
efficient or difficult to incorporate. The resulting images, 
formed by inverting reciprocal space diffraction patterns, con¬ 
tain quantitative information that encodes local physical pa¬ 
rameters such as permittivity, density, and atomic displace¬ 
ment at sub-beam-size spatial resolutions. 

When implemented with hard x-rays, these coherent 
diffraction imaging (CDI) techniques have enhanced our un¬ 
derstanding of the internal structure of nano- and meso-scale 
materials, especially in operating environments that are dif¬ 
ficult to access with other probes. Furthermore, x-ray mi¬ 
croscopy methods based on Bragg diffraction are of particu¬ 
lar interest because the sensitivity of x-rays to crystalline dis¬ 
tortions in materials can be leveraged to reveal the interplay 
between structure and properties without disturbing environ¬ 
mental boundary conditions. However, the routine application 
of inversion methods to coherent hard x-ray Bragg diffrac¬ 
tion is still limited by stringent experimental requirements and 
long measurement times. 

Given the potential impact of non-destructive 3D structural 
microscopy and the limitations of current 3D Bragg coher¬ 
ent x-ray inversion methods, advances in Bragg phase re¬ 
trieval methods that facilitate the rapid imaging of crystal lat¬ 
tice behavior in realistic environments are critically important. 
Here, we introduce a new coherent Bragg diffraction imaging 
approach, three dimensional Bragg projection ptychography 
(3DBPP), that provides such a capability. 3DBPP enables 3D 
image reconstruction from a series of 2D Bragg diffraction 
patterns measured at a single incident beam angle, thus form¬ 
ing a new mode of inversion-based 3D strain-sensitive imag¬ 
ing. As we discuss in this article, 3DBPP is a hybrid real / 
reciprocal space technique that takes advantage of the high 
angle of separation between the incident and diffracted beam 


in a Bragg diffraction geometry and eliminates the need to 
change the sample angle. In this way, 3DBPP fully exploits 
the 3D information encoded in a set of 2D Bragg diffraction 
patterns and paves the way for new high-throughput in-situ 
nanomaterials imaging studies. 


BACKGROUND 

In this work, we demonstrate a method by which fixed- 
angle 2D coherent diffraction patterns from a focused-beam 
scanning probe measurement are used to image 3D lattice dis¬ 
tortions in nanoscale crystalline structures. One key develop¬ 
ment that we present here is a new reconstruction algorithm 
that enables this new 3D sttTictural imaging capability. This 
algorithm features concepts utilized in other reconstruction- 
based microscopy methods, namely ptychography and tomog¬ 
raphy. In order to put our 3DBPP method in context, we 
briefly cover the essential concepts of Bragg coherent diffrac¬ 
tion imaging, ptychography, and tomography that relate to our 
new approach. 

X-ray microscopy of nano- and microscale crystals using 
coherent Bragg diffraction imaging allows three dimensional 
internal strain fields to be quantitatively measured GlEl. In its 
first implementation, Bragg coherent diffraction imaging in¬ 
volved illuminating a sub-micron-sized isolated crystal with a 
much larger coherent x-ray beam and measuring the diffracted 
intensity with an area detector Q. The far-field intensity 
distribution in the vicinity of a Bragg peak was recorded by 
scanning the fully illuminated crystal through the Bragg rock¬ 
ing curve (changing the sample angle relative to the incoming 
beam), combining successive 2D slices into a 3D Bragg inten¬ 
sity pattern. A three dimensional image of the diffracting crys¬ 
tal can then be numerically retrieved from the measured inten¬ 
sity distribution using iterative algorithms that constrain the 
extent of the crystal to within a support that roughly matches 
its size and shape ilia. This phase retrieval approach has 
been successfully used for in-situ 013 and ultrafast il imag- 
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ing studies of isolated micro- and nanocrystals. However, the 
intrinsic limits of the method |l2| constrain its application to 
a limited class of materials. For example, extended crystals 
or samples with large internal strain fields cannot readily be 
imaged with this approach. 

Recently, x-ray Bragg ptychography methods have been de¬ 
veloped that eliminate the requirement for isolated crystals 
and accommodate a broader range of samples. Originally 
proposed for electron microscopy Ham and developed ex¬ 
tensively with x-rays in the transmission geometry iniiii, 
the present form of ptychography consists of inverting a set of 
far-held diffraction intensity patterns collected from overlap¬ 
ping regions of the sample illuminated with a localized beam. 
Thus, specihc regions of interest can be imaged in continuous 
samples. As in the isolated crystal method, 3D images have 
been reconstructed from a 3D Fourier space intensity distri¬ 
bution about a Bragg peak nma. Although this approach 
provides access to the internal strain distribution in targeted 
regions of nano-structured crystals Ha, its implementation 
requires the measurement of diffraction patterns over a full 
angular rocking curve at each overlapping scan position, re¬ 
sulting in long measurement times and uncertainties with re¬ 
gard to registration. 

In general, the structural information encoded in a single 
far-held coherent diffraction area detector measurement is re¬ 
lated to the 3D structure of the sample. In the far-held regime, 
the components of the diffracted held exhibit a 3D depen¬ 
dence that is related to the illuminated scattering volume by 
a Fourier transform. An x-ray area detector accesses a two- 
dimensional “slice” through the intensity of this 3D recipro¬ 
cal space distribution. The slice, dehned by the plane of the 
detector, corresponds to the squared magnitude of the Fourier 
transform of a set of line integrals through the scattering vol¬ 
ume along the direction of the exit beam wave vector. In other 
words, the slice is the Fourier transform of a 2D projection 
of the illuminated sample volume. This property of x-ray 
diffraction is known as the “slice-projection theorem” (SPT) 
M Sec. 6.3.3], and it enables 2D and 3D imaging by various 
methods in both transmission and Bragg reflection geometries 

CMISl. 

In particular, the slice-projection theorem forms the foun¬ 
dation of modern x-ray computed tomography (CT) imaging 
algorithms. At its core, CT involves measuring projections of 
the scattering volume as a function of azimuthal sample angle. 
Whether obtained by direct or inversion-based methods, each 
real-space projection is Altered, then propagated along the di¬ 
rection of integration, a process known as backprojection. The 
Altered backprojections at each angle, which represent each 
projection in three dimensions in a manner consistent with the 
SPT, are summed to yield a 3D image E^ . However, the 
quality of a CT reconstruction depends strongly on angular 
diversity, and projections must be recorded over nearly the en¬ 
tire range of sample orientation 120°) to achieve a robust 
tomographic reconstruction imiini, iMl Sec. 6.2]. 

Below, we introduce a new inversion-based 3D structural 
imaging method that retrieves the sample without requiring 


angular diversity. This method relies on the iterative mini¬ 
mization of a cost-function that relates a set of measured 2D 
high-angle Bragg coherent diffraction patterns to a three di¬ 
mensional real-space sample structure. 

PRINCIPLES OF THE 3D BRAGG PROJECTION 
PTYCHOGRAPHY 

The key insight of this paper is that, by virtue of the ge¬ 
ometry required to satisfy a crystalline Bragg condition, the 
angular diversity critical for transmission-geometry CT can 
be replaced by translational diversity of a localized beam. We 
demonstrate that a 3D image can be obtained from only one 
slice of the rocking curve by presenting numerical and ex¬ 
perimental demonstrations as well as an analytical derivation 
found in the Supplemental. 

At a Bragg condition, the incident and scattered beam direc¬ 
tions are not collinear, but in the hard x-ray regime, they are 
typically separated by tens of degrees. As a result, the spatial 
information that is collapsed along the exit beam direction in 
any individual 2D projection can be encoded by translating a 
localized x-ray beam in the plane normal to ki, the propaga¬ 
tion direction of the incident beam (Figure a, b)). We use 
this principle as the basis of our new three dimensional Bragg 
projection ptychography technique, and we demonstrate that 
the 3D structure of a diffracting crystal is encoded in, and 
can be reconstructed from, a set of fixed-angle 2D coherent 
Bragg peak intensity patterns using spatial diversity. These 
patterns are collected by scanning the position of a localized 
x-ray beam relative to the crystal such that the beam footprint 
overlaps between neighboring beam positions. Thus, spatial 
oversampling is introduced, enabling the use of ptychographic 
phase retrieval methods while ensuring sensitivity to compo¬ 
nents of the 3D structure due to the Bragg geometry. 

In 3DBPP, the phases of the diffracted held are retrieved by 
defining a cost function and gradient that relate the 3D struc¬ 
ture of a crystal to the observed 2D Bragg diffraction inten¬ 
sity patterns. At a given beam position, the intensity pattern 
recorded by the x-ray camera in the far held is the squared 
amplitude of the diffracted wave held j: 

(Ij) = ( 1 ) 

where Ij is the intensity pattern corresponding to the j-th 
beam position. The intensity pattern recorded in the detec¬ 
tor is subject to counting statistics, and so Equation [T] defines 
an expectation value {Ij). The diffracted wave T'j is related 
to the 3D diffracting crystal p and the beam Pj (also referred 
to as the probe) according to Il25ll26l : 

T-, = PUP.p, (2) 

where P is the (bidimensional) Fourier transform, IZ is the 
(3D) x-ray projection operator along the exit beam direction 
ll24l Sec. II.2]. In this Equation, Pj and p are complex¬ 
valued three dimensional quantities. In the Bragg geometry. 
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a) Cross section of numerical sample p 



Internal density voids 





FIG. 1. The principles of 3DBPP. a) To illustrate the different 3DBPP reconstruction operators, a 3D strain free crystal was generated 
(p) containing two internal voids, and a 3DBPP experimental geometry was simulated (b) using a focussed gaussian beam profile at a high 
angle Bragg condition defined by the ki and kf incident and exit beam vectors. We define real space axes {rx,ry,rz) such that kr || rz and 
{rx, Ty) lie in the detector plane, conjugate to {qx, qy) in the far field. In order to calculate diffraction patterns, projections along kf of the 
illuminated crystal (TZPjp) in the {rx,ry) plane are determined (c) at beam positions that intersect the two voids in p. The far-field diffraction 
amplitudes IfiTjl determined by Fourier transforming the projections are shown in d). The inverse process is shown in e) where, starting from 
the far-field diffraction 'Tj , the inverse Fourier transform again yields the 2D projections of the illuminated crystal. The operator critical to 
3DBPP is the back-projection (7?.^) of the quantity along the kf direction. The backprojection “stretches” the 2D projections along 

kf within a 3D support. Here, the support is made of a pair of planes that limit the extent of the reconstruction in the 2 direction. Here, the 
backprojection intersects the probe at an oblique angle, effectively localizing the scattering volume. With 3DBPP, p was reconstructed (fig) 
using these operators from a set of 2D coherent Bragg diffraction intensity patterns (See Methods for details). 


the phase and amplitude of p are related to the Bragg struc¬ 
ture factor ini, which is sensitive to atomic-scale structure 
in the material. Our aim is to reconstract p, thereby imag¬ 
ing the diffraction stmcture factor in order to quantify local 
distortions of the crystal lattice. 

For 3DBPP phase retrieval, an inversion algorithm is used 
to reconstruct p from a set of J intensity patterns 
The numerical approaches that have successfully been used 
in transmission geometry ptychography experiments (i.e., dif¬ 
ference MAP Il27ll28l . ordered-subset (OS) / Ptychographic 
Iterative Engine algorithm (PIE) or gradient-type it¬ 

erations ll3TI43^ ) can be adapted to accommodate 3DBPP by 


incorporating a new gradient based on the cost function Q{p): 

j 

Q{p) = '^Qjip) with Qj(p) := II |«'j(p)| - V^lp. 

i=i 

( 3 ) 

In 3DBPP, this cost function yields the following gradient dj 
for each probe position; 

J = (4) 

where is the conjugate operator and TZ^ is the adjoint (i.e. 
















4 


backprojection) operator ll24l Eq. 2.31] associated with the 
forward-projection operator TZ. In this work, we incorporated 
the gradient in Equation]^ into an iterative OS/PIE algorithm 
that progressively lowered the cost function Q and recon¬ 
structed an image of the diffracting crystal p. (See Methods 
for more details.) 

As opposed to previous 3D Bragg ptychography ap¬ 
proaches iniiia, this new method is not one of pure Eourier 
synthesis in which the properties of the 3D reconstruction de¬ 
pend strictly on the measurement and sampling of a 3D vol¬ 
ume of reciprocal space. Rather, the 2D structural informa¬ 
tion of the sample that is encoded in the q^, Qy detector plane 
is used together with the real-space scanning probe positions 
to reconstruct the 3D sample structure. Thus, a successful 
3DBPP reconstruction must simultaneously retrieve two re¬ 
lated spatial components encoded in the phases of the diffrac¬ 
tion patterns. Eirst, each individual intensity pattern, when 
phased, gives access to the projected scattering volume for 
each beam position (the quantity in Eigure 1). Sec¬ 

ond, the relative phase relationship between the different in¬ 
tensity patterns governs the intersection of each backprojec¬ 
tion with the probe at each measurement position. This inter¬ 
section contains the spatial information along the projection 
direction that cannot be retrieved from an individual pattern. 
To aid in determining the latter phase relationship, various 
position-referencing constraints may be implemented depend¬ 
ing on the geometry of the sample. The thin film reconstruc¬ 
tions featured in this work were constrained with a support 
that confined the extent of the object along the surface-normal 
direction (as shown in Eigure 1). Eor other samples, differ¬ 
ent position-referencing constraints may be implemented that 
serve the same function. 

A numerical demonstration of 3DBPP is featured in Eig- 
ure We define a 3D thin-film object of uniform density, 
and introduced two cubic voids in the center of the film. A 
numerical set of 2D coherent diffraction patterns was gener¬ 
ated using Equation|^by scanning a purely real gaussian beam 
through the central region of the film in overlapping steps. 
Using 3DBPP implemented with an OS/PIE inversion algo¬ 
rithm, the 3D internal structure of the film in the scanned field 
of view was successfully reconstructed from this set of 2D 
patterns (Eigure [^f, g)), thus demonstrating the viability of 
our approach (see Methods for more details on this numerical 
test). 

EXPERIMENTAL DEMONSTRATION 

To demonstrate the efficacy and potential of 3D Bragg pro¬ 
jection ptychography, we imaged the internal lattice displace¬ 
ment field within a sub-micron-scale crystalline component 
of a semiconductor prototype device (Eigure |^. The device 
structure consisted of periodic embedded SiGe (eSiGe) crys¬ 
tals 460 nm in width that were epitaxially grown to a thickness 
of 65 nm interstitially between 60-nm-wide linear silicon-on- 
insulator (SOI) channels (shown in Eigure |^a)). This sam¬ 


ple design and processing (see Methods) resulted in a com¬ 
plex internal stress state and an accompanying strain field that 
evolves within the eSiGe stressors in the x and z directions, 
and is self-similar along y ll35l . In particular, the variation 
of the eSiGe strain field as a function of depth - variations to 
which 2D BPP is insensitive ll35]l - has a direct impact on the 
electron mobility within the SOI channels in this system and 
is a critical parameter for nanoscale device engineering. The 
3D reconstruction of an eSiGe strain field presented here via 
3DBPP demonstrates the power of this method for in-situ non¬ 
destructive structural imaging of functional crystalline nano¬ 
materials. 

Coherent nanodiffraction patterns were measured from sev¬ 
eral adjacent eSiGe stressor crystals with a zone-plate-focused 
hard x-ray beam at a symmetric 004 Bragg condition in which 
the angle between ki and kf was about 60° (see Methods and 
Eigure]^. We note here that the experimental data collection 
methodology of 3DBPP is exactly the same as that of pre¬ 
viously reported 2D BPP experiments ifThUTSl [35l . however 
using our new imaging concept, 3D images can now be re¬ 
constructed from data that previously yielded 2D projection 
images. Using these Bragg diffraction patterns, a 3DBPP re¬ 
construction was generated using an iterative algorithm that 
incorporated the gradient in Equation into an ordered sub¬ 
set / PIE framework lEolIiI]. Because the coherent diffrac¬ 
tion observed in the vicinity of the eSiGe 004 Bragg condition 
(denoted by the reciprocal space vector G 004 ) was well sepa¬ 
rated from scattering from the other components of the device 
(i.e. SOI and substrate), the reconstruction represents only the 
stressor structures. To aid in the determination of the relative 
phase relationship of the intensity patterns, a 90-nm-thick sup¬ 
port that confined the extent of the reconstruction along z was 
incorporated into the reconstruction. This support constrained 
the intersection of each backprojection and respective probe 
function to within a physically realistic pair of parallel planes 
surrounding the 65-nm-thick film (see Methods). 

The resulting eSiGe reconstruction is shown in Eigures 
Two linear SOI channels that are each 60 nm wide run along 
the y direction. The footprints of these channels are depicted 
in Eigure 3(c,d) by two parallel green strips. The SiGe recon¬ 
struction in our 3D field of view therefore features four verti¬ 
cal SiGe/SOI interfaces (along the edges of the green strips) as 
well as the interface beneath the SiGe stressor material. The 
outer boundary of the isosurface shown in Eigure 3(a) was de¬ 
fined by the field of view of our 2D beam position scan in the 
X, y plane. 

As is expected for strained crystals, the reconstructed 3D 
quantity p is a complex-valued function, p = \p\ eyzpi(j). The 
amplitude |p| is directly related to the electron density of the 
crystal and the phase (p is sensitive to a component of the dis¬ 
placement field in the crystal given by ^ = G004 • u, where 
u is the crystalline displacement field relative to an arbitrary 
crystalline reference 0. 

The external shape of the reconstructed amplitude is shown 
in Eigure as an iso-surface corresponding to 27% of the 
maximum amplitude value. The reconstructed shape shows 
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FIG. 2. Experimental geometry, (a) A schematic of the device architecture is shown featuring the embedded SiGe stressor geometry. For 
3DBPP, coherent diffraction patterns were measured at the specular 004 SiGe Bragg peak at a set of beam positions that effectively overlapped 
in a series of spiral patterns. The diffraction patterns from all 707 positions were phased simultaneously with 3DBPP. (b,c) Prior to diffracting 
from the eSiGe, the beam wavefront was reconstructed using Fresnel ptychography of a test pattern in the transmission geometry, (d-f) 
Examples of 004 coherent Bragg intensity patterns used for 3DBPP imaging are shown from different regions of the stressor indicated in (a). 
Patterns (d-f) correspond to the positions indicated with blue dots in (a) from left to right. See Methods and Ref for more details. 


the eSiGe features expected in the scanned field of view, in¬ 
cluding the two missing strips in amplitude corresponding to 
the positions of SOI channels. The internal structure of the 
reconstructed stressors is shown in Figure [^c,d). The verti¬ 
cal cuts through the reconstruction show that the amplitude, 
which is closely related to the crystal density, is mostly ho¬ 
mogeneous with well defined in-plane edges where the eSiGe 
meets the SOI (further discussion below). However, in this 
case, the unique strength of Bragg peak inversion lies in visu¬ 
alizing local lattice distortions in a crystal by way of the recon¬ 
structed phase. For this stressor, the displacement of the crys¬ 
tal lattice along the [001] direction, Uqoi. was derived from 
the phase of the 3DBPP reconstruction (uooi = ^i>/|Goo 4 |), 
and is shown in Figure [^d). As expected for this device ar¬ 
chitecture, Uqoi varies in the (x, z) plane. The key advantage 
of 3DBPP is demonstrated by the fact that out-of-plane lattice 
distortions in this plane that require a 3D reconstruction were 
resolved, and that this was accomplished with a fixed-angle 
scanning probe diffraction measurement. 

A crucial issue in this study is the quantification of the res¬ 
olution obtained in the reconstruction, as this will dictate the 
applicability of 3DBPP towards the measurement of lattice 
distortions in different nanoscale crystalline systems. To ac¬ 
complish this, we used a method similar to that used in Ref 
ll^ . which is based on an analysis of the complex recon¬ 
structed density p in terms of spatial frequencies that persist 
above noise, adapted for a continuous extended 3D sample 
(see Methods). In the plane of the detector, anisotropic spa¬ 
tial resolutions of ^18 nm along x and ^17 nm in the detec¬ 
tor direction orthogonal to x were obtained, while along the 
projection direction (i.e. along kf), the resolution was ~25 
nm. The differences in these resolutions serve to highlight 


the resolution anisotropy that can occur with this technique, 
though we note that sub-beam-size resolution were achieved 
in all three directions (see Supplemental). As in other pty¬ 
chography experiments, better knowledge of the beam profile 
and beam positions yields a more accurate and higher reso¬ 
lution 3DBPP image. Strategies to mitigate uncertainties in 
these quantities can potentially be adapted from transmission 
geometry ptychography to 3DBPP in order to improve image 
fidelity. 

The structural fidelity of the method was evaluated by com¬ 
paring a cross-section of the reconstruction with a linear elas¬ 
tic boundary element method (BEM) model of the eSiGe 
stressor under the nominal mechanical boundary conditions 
of the device IIJTI . Figure [^(a-c) shows the amplitude, phase, 
and displacement field maps from an (x, z) cross-section of 
the central reconstructed stressor (cut-plane A in Figure 
alongside those of the corresponding BEM model. The BEM 
model shows that the out-of-plane component of the displace¬ 
ment field Uqoi evolves within the (x, z) plane, intensifying 
at the top corners to a value of approximately —2.5 A relative 
to the center (Eigure|^(f)). This variation in Uooi is the result 
of an elastic response of the SiGe lattice due to a change from 
a near-biaxially stressed state at the center of the stressor to 
a more uniaxially stressed state at the eSiGe / SOI interface. 
The model Uooi displacement was converted to a complex 
crystal density pbem calculated for the 004 Bragg diffrac¬ 
tion condition, and a Eourier filter was applied to determine 
phase and amplitude distributions {Pbem^ These model dis¬ 
tributions are consistent with the experimental Bragg diffrac¬ 
tion signal levels (see Methods). As seen in Eigurej^d), the 
Eourier filter modified the rectangular cross-section of the am¬ 
plitude because high-spatial-frequency components not exper- 
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FIG. 3. 3D Bragg projection ptychography experimental results Coherent diffraction patterns measured at the 004 SiGe Bragg peak were 
phased using 300 iterations of the OS/PIE algorithm adapted for 3DBPP. (a) The cost function, Q{p), is shown as a function of iteration 
number, (h) The isosurface of the amplitude of the resulting reconstruction, showing gaps at the positions of the SOI channels. The internal 
structure of the SiGe material in the field of view is revealed by way of cuts through the isosurface depicting the amplitude of the reconstruction, 
which is closely related to the material density (c), as well as the crystal lattice displacement along the [001] direction (d) which is derived 
from the phase of the reconstruction. 


imentally observed in our measurement were filtered out ED. 
In the center of the amplitude remains homogeneous, 

but the amplitude envelope falls off near the top comers where 
the internal displacement field of the crystal is expected to 
vary more rapidly. As a result of this displacement field, 
the phase evolves by ~ 27r in the top corners of the stres¬ 
sor before the amplitude envelope drops off. These features in 
Pbem’ which come about from the predicted internal lattice 
behavior of the stressor as well as the limited dynamical range 
in the measured intensity, are all observed in the experimen¬ 
tal 3DBPP reconstmction (figure [^,b). This imaging experi¬ 
ment establishes the effectiveness of our method in using 2D 
diffraction patterns to reconstract relevant 3D structural infor¬ 
mation in targeted regions of nanoscale crystals - a capability 
that can be used to explore and discover new phenomena over 


a broad range of complex strained crystals. 


DISCUSSION & CONCLUSION 

Our experimental results demonstrate that 3D images of 
strained crystals can be reconstructed without angular diver¬ 
sity at resolutions appropriate for nanoscale imaging. How¬ 
ever, certain factors must be considered for a successful exper¬ 
iment. Chiefly, the effectiveness of localizing sample structure 
along the beam propagation direction (via the Pj and TZ^ oper¬ 
ators) diminishes as the angle between ki and kf approaches 
zero, while it is maximized at 90°. Secondly, the image res¬ 
olution along the backprojection direction is coupled to both 
the beam size and the step size, while the resolutions parallel 
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FIG. 4. Comparison with linear elastic model: A cut through the eSiGe stressor reconstruction is shown as both amplitude (a) and phase 
(b), with the edges of the 90-nm-thick support planes indicated. The cross section shown corresponds to cut-plane A in Figure [^c). In (c), 
the phases in (b) were converted to units of lattice displacement along the [001] direction (uooi) confined within an envelope representing the 
nominal edges of the 65-nm-thick film. For comparison, a boundary element method (BEM) linear elastic model of uooi in the stressor (f) 
was converted to the complex density expected for 004 Bragg diffraction, pbem- The amplitude and phase of pbem and p^g'j^j^f are shown in 
(a,b) before and after application of a low-bandpass Fourier filter that replicates the experimental signal level. 


to the detector plane are limited by the extent of the reciprocal 
space signal. We note that because the method involves suc¬ 
cessive backprojections, the reconstruction is susceptible to 
noise amplification that degrades the image at high iteration 
numbers in a manner analogous to CT reconstructions with 
noisy data (see Supplemental). 

Efficient three dimensional microscopy of strain fields in 
targeted regions of nanoscale crystals under balanced, undis¬ 
turbed boundary conditions will enable powerful new stud¬ 
ies of in-situ materials behavior. By utilizing a new inver¬ 
sion strategy that efficiently uses Bragg diffracted photons 
from a simpler and dose-efficient experiment, three dimen¬ 
sional Bragg projection ptychography is a significant step in 
the development of 3D x-ray microscopy, especially with re¬ 
gard to radiation-sensitive materials. As compared to Bragg 
ptychography with rocking curves, 3D images of crystals can 
be reconstructed using a single diffraction pattern measure¬ 
ment at each point, providing a 50- to 100-fold reduction in 
dose and acquisition time. Furthermore, our approach simpli¬ 
fies the experimental requirements such that 3DBPP can be 
implemented on a wide range of coherent synchrotron x-ray 
beamlines. More broadly, the integration of projection and 
backprojection operations into a phase retrieval algorithm can 
potentially be applied to address challenges in transmission 
geometry computed tomography, including the alignment of 
radial projections with incoherent CT and potentially enabling 
global 3D inversion of tomographic far-field coherent diffrac¬ 
tion patterns. 


METHODS 

Numerical 3DBPP demonstration: A numerical sample 
in the shape of a rectangular cuboid (Figure 1) was generated 


that was 40 x 30 x 190 pixels in dimension with two cubic 
voids in the sample with respective edge lengths of 7 and 17 
pixels. The numerical object (denoted as p) was purely real 
and had a uniform density outside of the voids. Figure 1(a) 
shows a cross-section through the middle of p, showing the 
two voids in density. The isosurfaces in Figure l(b,c,f) all 
depict the density from one outside edge of the object to the 
central cross-sectional plane shown in (a). Though the entire 
numerical object is not depicted in the isosurfaces, this visual¬ 
ization helps emphasize and clarify the position and shape of 
the voids before and after 3DBPP reconstruction. The object 
is symmetric about the plane in (a). 

The numerical probe used for this example was purely real 
with a gaussian amplitude profile with a full width at half 
max width of 4 pixels. To simulate a high-angle Bragg condi¬ 
tion similar to the eSiGe experiment, the incident angle of the 
probe and the angle of the exit beam (kj and kf) were both 
set to 30° with respect to the top surface of the object. In or¬ 
der to generate a simulated data set, the position of the probe 
was scanned through the object in a 5 x 61 point grid with a 
two pixel step size. At each point j, coherent Bragg intensity 
patterns were calculated according to (Ij) = The 

resulting noise-free intensity patterns were phased using the 
OS/PIE 3DBPP algorithm described in the text to generate a 
reconstruction. The support used in the reconstruction (yellow 
surfaces in Figure 1) consisted of parallel planes separated by 
40 pixels, matching the dimension of the original object in 
the z direction. This support acted only to constrain the top 
and bottom surfaces of the object, as in the experimental eS¬ 
iGe reconstruction. The final reconstruction in Figure 1(g) 
was generated after 200 iterations of the 3DBPP algorithm. 
The resolution of this reconstruction along kf was found to be 
2.6 pixels, demonstrating sub-beam-size resolution of 3DBPP 
along the projection direction (see Supplemental). 
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Sample preparation: The sample was a lithographically 
prepared prototype device lIMl . consisting of a series of 460- 
nm wide trenches etched into a silicon-on-insulator (SOI) 
wafer. The trenches were subsequently filled with epitaxial 
embedded Sio. 8 Geo .2 (eSiGe) stressors that were 65 nm thick, 
as confirmed by transmission electron microscopy. The ge¬ 
ometry and mechanical boundary conditions of these epitax¬ 
ial device components are such that strain gradients are ex¬ 
pected to build up near the vertical interface between the eS¬ 
iGe and SOI channel regions. The device SOI/eSiGE archi¬ 
tecture shown in Figure]^ a) repeats for > 100 periods in the 
X direction, and is invariant and self-similar in the y direction 
for tens of microns. 

Data collection: The coherent diffraction measurements 
were done at the Hard X-ray Nanoprobe beamline IllHIl. 
A Fresnel zone plate was used to focus 9 keV energy x-rays 
with wavelength A = 0.137 nm to a ~40 nm diameter focal 
spot (FWHM intensity of the central focus peak). Prior to col¬ 
lecting nanodiffraction patterns from the eSiGe stressors, the 
complex wavefront of the beam was determined in the focal 
plane with transmission-geometry Fresnel ptychography us¬ 
ing a test object ll42l |4^ . (Figure |^b-c)). Using wavefront 
propagation, the wave held of the beam about the focus was 
calculated in order to determine the 3D probe function in the 
Bragg geometry. 

For 3DBPP imaging, coherent nanodiffraction patterns 
were measured with an area detector at the 004 eSiGe Bragg 
condition with 10 second exposures. In this symmetric 
diffraction geometry the incident and exit wave vectors, ki 
and kf, were separated by 59.5°. The beam was scanned nor¬ 
mal to ki in a series of spiral-like patterns with an effective 
step size of ~13 nm separating adjacent points, dehning an 
effective held of view of ~ 900 x 200 nm on the surface 
of the device. In total, 707 coherent Bragg diffraction pat¬ 
terns were phased, covering a sample area shown in Figure 
l^a). Examples of coherent Bragg diffraction patterns from 
this data set measured at different positions of a single stres¬ 
sor are shown in Figure|^d-f). Further details can be found in 
Ref 051 . especially with regard to our implementation of the 
spiral patterns employed in order to mitigate beam damage. 

Inversion procedure: The reconstruction was initialized 
with a 3D starting guess of p consisting of random real values. 
The beam wavefront reconstructed from test pattern data (Fig¬ 
ure l^b)) was used to generate a 3D incident focused-beam 
probe P that was kept constant over the course of the recon¬ 
struction. A total of 300 iterations were performed for the 
3DBPP phase retrieval. A support consisting of two parallel 
planes separated by 90 nm was applied to each new update of 
p such that amplitude outside the volume between the planes 
was set to zero. A /3 value of 0.8 was used during the OS/PIE 
iterations OOl . Similar to the approach in Reference 1311 . we 
utilized a pre-conditioning matrix to aid in the convergence of 
the inversion. 

Estimation of the spatial resolution. Following an idea 
proposed in Reference ll^ . the experimental SiGe recon¬ 
struction (p) is used to calculate a set of diffraction intensity 


patterns using the known probe and probe positions. In these 
diffraction patterns, intensity below a threshold level of T ~ 1 
photon was set to zero resulting in a set of thresholded inten¬ 
sity patterns The spatial resolution of 'p is esti¬ 
mated using by minimizing the following monodi¬ 

mensional criterion, which is closely related to the criterion 
expressed in Equation but which uses the thresholded in¬ 
tensity patterns rather than experimentally observed ones: 


HII “ l-^^^jPa(q)l IT- 

i=i 

In this expression, 'p^ is a filtered (low-resolution) version of 
the reconstruction 'p defined by its 3D Fourier transform: 


{J^3DPa) (q) 


0 if |(.F3r,;0a)(q)| < a 

{J^ 3 Dp){q) otherwise. 


The a-hltered version of the reconstruction yields a 3D image 
that excludes some of the high spatial frequencies. This crite¬ 
rion is minimized at a value of a = acr, and acr is then used 
to determine the spatial resolution of 'p. We take the Fourier 
intensity \p3D'PaaA^ determine the extent of the non-zero 
reciprocal space envelope (qmax) spanning from the central 
pixel along each of three directions: the two directions parallel 
to the plane of the detector and in the direction of projection. 
Finally, qmax in each of these directions was converted to an 
estimate of resolution using the relation: r = 2Tr/{2qmax)- 


Fourier filtered BEM model: At an exposure time of 10 
seconds, the experimental Bragg coherent diffraction signal 
from the eSiGe stressors fell below the noise level of the 
detector at a signal level of ~0.25% of the maximum ob¬ 
served intensity, below which diffracted signal was not de¬ 
tected. This effectively removed high-frequency components 
of the diffracted Bragg peak measurement at each position. 
Because of the hnite measured signal, our experimental re¬ 
construction can be thought of as a low-bandpass filtered ver¬ 
sion of the reconstruction that would have been obtained with 
noise-free high-dynamic-range data. 


In order to compare the results of the linear elastic model 
with our experimental reconstruction, a comparable low- 
bandpass hlter was applied. The BEM model, when converted 
to pbem for the 004 Bragg condition (Figure 4(d,e)), contains 
high spatial frequency features (sharp edges, rapidly varying 
phase ramps) that would not all be detected at the signal level 
corresponding to our experiment. To represent pbem in a 
manner that was consistent with the hnite experimental signal, 
frequency components of the Fourier transform of pB em that 
fell below 0.25% of the maximum intensity of the FT were 
zeroed. The inverse Fourier transform of this quantity yielded 
Pbem’ which is an estimate of a 3DBPP reconstruction of the 
BEM model that is consistent with the signal-to-noise ratio of 
our experiment. 
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Supplemental Section: 

Derivation of 3D structural encoding in 3DBPP 

As introduce above, the electronic density of the diffract¬ 
ing crystal and the probe function are denoted, respectively, 
p{r) and P{r). For sake of simplicity, we consider the bidi- 
mensional case depicted in Figure 1(a). This simplihed con¬ 
struction more clearly conveys the fact that structural informa¬ 
tion in the direction is encoded in a scanning probe Bragg 
diffraction experiment without a rocking curve. Our conclu¬ 
sions remain unchanged when the system is extended to 3D. 

The orthonormal reference frame in the real space (e^,, e^) 
is chosen so that it is the conjugate of the k^) reference 
frame, and so that k^ is aligned with the detector. By defini¬ 
tion, any r (in real space) and q (in reciprocal space) reads 

r = and q = q^k^ + q^k^, 

with ||ea;^ 2 || = ||fca;,z|| = 1- 

Let us consider the displacement of the probe along the di¬ 
rection Ez in discrete steps: 

-Pm(f) := P{'i' — rnAez), where m = 0, ■ ■ ■ M — 1. 

Here, the beam step size is defined as A S R. In practice, the 
displacement of the probe will also have a component along 
e^, but we neglect this component here because it does con¬ 
tribute to the encoding of information along e^. 

As a result of the Slice Projection Theorem ll24li rSec. 


6.3.3], when a Bragg condition is satisfied, the coherent far- 
held intensity pattern observed in the detector is given by: 

= \'^miqz = 0,qx)\^ (5) 

where ■!/)„ is the Fourier transform (indicated by the “ ~ ” 
symbol) of the m-th sample exit-held: 

il^mir) := P{r - mAez)p{r) ( 6 ) 

For the purposes of this derivation, we assume that the phases 
of ipm are known {i.e. an appropriate phase retrieval strat¬ 
egy was implemented, such as the one we presented in the 
manuscript). Thus, we dehne the following quantity: 

Cim,qx) ■=i’miqz =0,qx)- (7) 

(^{m,qx) represents a set of complex-valued far-held coher¬ 
ent diffraction wave helds at M different probe positions ob¬ 
served at a hxed sample angle (no rocking curve). Similarly, 
in real space I’x) represents the corresponding set of exit 
wave functions at the sample at each probe position. The ques¬ 
tion we seek to answer is whether structural information about 
the sample along is encoded in either of these quantities. 

To answer this question, we begin by expressing Equation 
(|^ as: 

il^ 7 n{r) = p{r) X [P(r) (g) S{r - mAe^)], 

where g) is the convolution operator. Taking the Fourier trans¬ 
form leads to: 


= Kq) ® \Hq) x 


( 8 ) 


= g-127r(mA)(3j 



')P{q - dq' 


(9) 


We now assume that p and P are both separable functions, i.e., 


p{r) = Pz{rz) X px(rx) p{q) = Pz{qz) x Px{qx) 

P{r) = Pzivz) X Px{rx) ^ P{q) = Pz{qz) x Px{qx) 


( 10 ) 

( 11 ) 


This assumption allows us to most clearly demonstrate how spatial information along can be encoded without perform¬ 
ing a rocking curve. To do this, we express Equation (|^ as: 


= (^Px ® Px'^ {qx) X 


Pz{q'z)Pz{qz-q'z)eP^^^^^P dq'z 




( 12 ) 
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According to Equation Q, at a single Bragg angle we can 
only access the far-field diffracted wave field at = 0. The 
reciprocal space information we collect as a function of probe 
position can therefore be expressed as: 

C(to, Qoo) = {px ® Px^ (Qx) X hx{rx = mis) (13) 

were hz is the inverse Fourier transform of hz (q) = pz (<?) x 
Pz{—q). Another convenient expression for Equation ([T^ is: 

C(m, qx) = [px ® Px'j (Qx) X {pz ® P'z) {tz = mis) (14) 

where P'zivz) := E’z(—^z). Finally, we can derive the depen¬ 
dence in real space as a set of sample exit waves by taking the 
inverse Fourier transform along q^'- 

C(m,rx) = Px(rx) x Px(rx) x (pz ® P'z) (rz = mA). (15) 

From this expression, we see that the sample exit wave held, 
when measured as a function of probe position m, depends on 
Tz via the convolution of pz with P),. 

Hence, the structural information contained along the rock¬ 
ing curve direction is encoded in a set of single-angle diffrac¬ 
tion patterns when they are measured with a scanned focused 
beam. This conclusion also holds true when a third dimen¬ 
sion (into the page) is added. This is because the structural 
information about the sample along Vy is directly encoded in 
reciprocal space with an area detector that can access q^ and 
qy, and does not change the discussion about sample structure 
along Tz. 

Comparison of Bragg and transmission geometries 

In order to emphasize the role of the Bragg geometry in 
3DBPP, we present a comparison between images recon¬ 
structed from diffraction patterns simulated in both the Bragg 
and forward scattering geometries (Figure SI). In both cases, 
intensity patterns give access to a 2D projection of the illu¬ 
minated volume. However, the relationship between the exit 
beam vector kf (the projection direction) and the incident 
beam direction ki in the two cases leads to very different three 
dimensional reconstructions. 

As outlined in the main text, in the Bragg geometry ki and 
kf are not colinear. Thus, scanning the localized incident 
beam in the plane normal to its propagation direction with 
overlapping steps enables the 3D structure of the sample to be 
reconstructed using the 3DBPP phasing approach. In this nu¬ 
merical example, the phases of the far field intensity patterns 
at each position are known a-priori. Therefore, a summation 
of the back-projections at each scan point multiplied by each 
respective probe can be done: PjP^j- This sum¬ 

mation of localized backprojections is shown in Figure Sl(b) 
for the numerical sample described in Figure 1 of the main 
text. The resulting image shows that the three dimensional 
sample structure can be recovered from a set of 2D projec¬ 


tions provided that ki and kf are separated by a non-zero 
angle typical of hard x-ray Bragg diffraction. This process 
yields an image that is a blurred version of the original ob¬ 
ject, and demonstrates how 3D structural sample information 
is encoded in a high-angle Bragg ptychography data set with¬ 
out requiring angular diversity. Starting with the same numer¬ 
ical diffraction information in the form of intensity patterns, 
3DBPP was used to recover the phases and reconstruct the 
object. As shown in Figure Sl(c), the 3DBPP reconstruction 
is sharper than the summation of localized back projections, 
and a sub-beam-size resolution is achieved (see next section 
in Supplemental) as has been shown in transmission-geometry 
ptychography experiments ll46l . 

As compared with the Bragg geometry, kj and kf are 
mostly colinear in a scanning beam transmission geometry 
diffraction experiment (the maximum angle between ki and 
kf is typically < 2°), as shown in Figure S 1(d). As a result, no 
structural information along the propagation direction of the 
incident beam is encoded in a 2D diffraction pattern measured 
in this geometry. This is demonstrated by the fact that both 
the summation of localized backprojections and the 3DBPP 
reconstruction of the numerical test object yield images that 
show the projected density profile of the object stretched along 
the kf direction. In other words, only the backprojection of 
the original object within the scanned held of view was re¬ 
constructed. In order to recover the three dimensional struc¬ 
ture of the sample from such a scattering geometry, scans must 
be done as a function of sample orientation and computed to¬ 
mography algorithms must be used for 3D imaging ifT^ . 

Three-dimensional dependence of resolution 

As in other coherent diffraction imaging techniques, the 
resolution of a 3DBPP reconstruction varies along different 
directions. In the plane perpendicular to kf, the resolution 
is primarily restricted by the extent of the diffraction enve¬ 
lope measured about the Bragg peak, which sets the effective 
numerical aperture in that plane. This type of g-range reso¬ 
lution limitation is typical of diffractive imaging microscopy 
techniques, and accurately describes the experimental resolu¬ 
tion of the eSiGe reconstruction along the axes of the detector 
plane. 

However, the factors that limit image resolution in the pro¬ 
jection direction are less obvious. The resolution in the kf di¬ 
rection does not stem simply from reciprocal space sampling, 
but is influenced by the scanning probe nature of the diffrac¬ 
tion experiment. In order to investigate the dependence of res¬ 
olution along kf as a function of scattering angle, beam size, 
and beam step size; 3DBPP reconstructions were done using 
calculated noise-free data from the numerical object shown in 
Figure S4(c). Various symmetric Bragg peaks were simulated 
with 29 = Zkfkf = 50,60,70,90°, and an exact support 
was used that conformed to the outer edges of the sample, 
leaving the voids to be reconstructed by the algorithm. Real 
gaussian beams with a FWHM of 4 and 6 pixels were used 
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FIG. 5. Suppl (SI). Bragg vs. transmission geometry 3DBPP (a) The Bragg geometry is simulated for 3DBPP of the test object featured in 
Figure 1 in the main text, (b) Because the phases of each diffraction pattern are known in this simulation the quantity j can 

be calculated, referred to as the summation of localized backprojections. (c,d) The result of 3DBPP phase retrieval using simulated noise-free 
diffraction patterns is shown, e) The same test object is scanned with the same focused beam in a transmission geometry, (f) The summation 
of localized backprojections does not accurately represent the voids, nor does the 3DBPP reconstruction (g,h) because 3D sample information 
is not encoded in a transmission ptychography experiment. 
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2 theta 2 theta 

FIG. 6. Suppl. (S2). Dependence of resolution: 3DBPP resolution along the kf direction is shown as a function of step size, scattering angle, 
and beam size from reconstructions of the numerical object shown in Figure S4(c) at a series of symmetric Bragg reflections. Noise-free data 
was used with a tight support that was the exact outer shape of the original object in order to isolate the effects of the other parameters indicated 
here. The solid lines are linear fits to each data set and serve to highlight the dependence of resolution on Bragg angle, step size, and beam 
size. 
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in this comparison. The resolution in the kf direction of the 
resulting reconstructions are shown in Figure S2. They are 
expressed as the FWHM of a gaussian blur filter kernel ap¬ 
plied to a set of amplitude line profiles along kf through the 
original object fitted to the corresponding amplitude profiles 
of the reconstruction. The mean value of a set of 15 line pro¬ 
file fits from each reconstruction are plotted with error bars 
of one standard deviation. Linear fits to the data are shown 
that are meant to be guides to the eye. As shown in the Fig¬ 
ure, sub-beam-size resolution can be achieved along the kf 
direction over an angular range of 20 ~ 40 — 90° for the con¬ 
ditions surveyed here. The resolution does, however, have a 
clear dependence on scattering angle, step size, and beam di¬ 
ameter. Though a description of the precise dependence and 
limits of 3DBPP resolution along kf requires further inves¬ 
tigation, based on this numerical study we conclude that the 
resolution of a 3DBPP imaging experiment is highest when 
the beam diameter and step size are minimized (though these 
parameters have limits in practice). Furthermore, the resolu¬ 
tion has a strong dependence on scattering angle, which we 
expect to be optimized when kf is normal to ki. As the 20 
angle approaches 0°, we encounter the situation described in 
Figure Sl(e-h), in which all structure along the kf direction is 
lost. 

In practice, the dependence and interconnectivity of the 3D 
resolution function of a 3DBPP reconstruction is undoubtedly 
a complex function of Bragg angle, beam size, beam ampli¬ 
tude and phase profile, beam step size, experimental position¬ 
ing accuracy and drift, signal-to-noise ratio, reconstruction al¬ 
gorithm, and other factors. Furthermore, the experimenter has 
considerable latitude in defining a support. In our experimen¬ 
tal reconstruction, the support was defined by a pair of parallel 
planes that were consistent with the expected thickness of the 
diffracting crystal in the out-of-plane direction. In the case of 
crystals that have thin-film-like geometries, we find that this is 
a suitable choice of support. In addition, the thickness of such 
crystals can be readily determined by other methods and can 
even be estimated from the diffraction patterns in the 3DBPP 
data set. We note that an effective support for a 3DBPP exper¬ 
iment primarily fulfills the role of a phase reference along kf, 
and thus can take on forms other than the ones we employed 
here. 

In Figure S3, we explore the dependence of 3D image fi¬ 
delity as a function of support size by way of numerical sim¬ 
ulation of a cubic crystal. Noise-free intensity patterns were 
generated from a 7 x 7 grid of beam positions illuminating a 
simple cube with an edge-length of 10 pixels. The angle be¬ 
tween ki and kf was 60°, as indicated in the Figure. The plot 
in Figure S2(a) shows the sum squared difference of pixels 
between the original cube porig and the resulting reconstruc¬ 
tion p determined according to \ \p — PorigW"^ as a function of 
support size. The support used in each case was cubic with 
edge lengths as shown in the plot. Two hundred OS/PIE it¬ 
erations were performed to generate a reconstruction, and re¬ 
constructions were done five times in each case. The mean 
I Ip — PorigW^ value is plotted with error bars of one standard 
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FIG. 7. Suppl (S3). Effect of support in 3DBPP. 3DBPP recon¬ 
structions of a cubic numerical test object were done using different 
sized cubic supports. At each support condition, five separate recon¬ 
structions were conducted. The mean final real-space object error 
Up — Porig IP is shown US a function of support size in (a), (b) Ex¬ 
amples of a reconstruction from each condition are shown (gray iso¬ 
surface) along with the corresponding support (yellow surface). The 
geometry of the experiment is conveyed by the incident and exit kj 
and kf vectors. 

deviation, and values of ||p — PorigW^ were determined ac¬ 
counting for the differences in origin of p relative to that of 
the Porig- As the edge length of the support increases from 
being exactly the size of Porig to being nearly twice its size. 
Up — PorigW^ increases and the reconstructions contain more 
artifacts that clearly affect the resolution. However, this test 
demonstrates that the cube reconstructions are recognizable 
even with a support that is nearly eight times the volume of 
the original object, suggesting that the choice of support is 
rather flexible in 3DBPP. 


Noise amplification during 3DBPP phasing 

In a transmission geometry, the numerical inversion of the 
x-ray projection operator TZ from experimental data, as im¬ 
plemented in CT for example, is highly sensitive to noise 
ll24l Sec. 4]. This sensitivity is a consequence of the SPT 
In the transmission geometry, the Fourier domain is “sliced” 
in polar coordinates, and high frequency components are sam¬ 
pled more sparsely as compared to low frequency components 
when the sample is rotated to measure the 3D intensity distri¬ 
bution. At a Bragg peak, on the other hand, the 3D intensity 
distribution can be measured with nearly parallel 2D slices as 
the sample is rocked in fine angular increments through the 
crystal rocking curve. This sampling of 3D reciprocal space 
at a Bragg peak is therefore done in a cartesian (though not 
orthogonal) frame lEl rather than a polar coordinate system. 
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As a result, the high and low frequency components of a 3D 
Bragg peak intensity distribution are sampled equally, and this 
leads to a robust ptychographic reconstruction when angular 
diversity is utilized in the Bragg geometry to measure the 3D 
intensity distribution at each probe position. 

For the 3DBPP technique introduced in this work, how¬ 
ever, it is not clear a priori if the method will be sensitive 
to noise because the 3D sample information is extracted via 
shifting overlapping probe positions rather than polar or carte¬ 
sian sampling of a 3D Fourier space far-held intensity distri¬ 
bution. Here, we addressed this issue by way of a numerical 
simulation. Coherent diffraction data was simulated from a 
3D numerical test object (Figure S4) similar to the one fea¬ 
tured in Figure 1 in the main text. While the ptychographic 
data was simulated in the same way as before, in this test the 
resulting intensity patterns were corrupted with Poisson noise 
in order to study its effects. From a starting guess of real ran¬ 
dom numbers, 100 iterations of OS/PIE were done in order 
to generate a starting estimate of the sample for the purposes 
of this test. From this starting estimate, 400 more iterations 
of the OS/PIE algorithm were performed, and the same ob¬ 
ject estimate was used as a starting point for 400 conjugate- 
gradient (CG) iterations. The CG inversion has been demon¬ 
strated to be convergent Sec. 5.2], and hence provides (in 
the asymptotic limit) a local minimizer of the cost-function 
Q given by Equation 4 in the main text. In both cases, the 
reciprocal space error metric Q was tracked as a function of 


iteration number as well as the real space object error metric 

Up ~ PorigiP- 

As shown in Eigures S3(a, b), both algorithms minimize 
Q, however their behavior is very different when evaluating 
how closely the reconstruction resembles the original object 
at each iteration in terms of ||p — PorigW^- At first, the CG 
algorithm improves the quality of the result more rapidly than 
the OS/PIE algorithm (Eigure S3(b)), however, after 52 it¬ 
erations the CG Up — PoriglP begins to increase, while Q 
continues to decrease. By contrast, with OS/PIE, the values 
of Up — PorigW^ and Q both decrease over the course of the 
400 iterations (however, this is not necessarily an indication 
of convergence). According to standard inversion theory ll45l 
Sec. 6], the behavior displayed by the CG algorithm suggests 
that 3DBPP is noise sensitive, meaning that high-resolution 
features in the object cannot be retrieved from the noisy data 
set without regularization. However, stable band-limited re¬ 
constructions are obtained in the early iterations by the CG, 
and similar behavior can be obtained by other gradient-based 
iterative algorithms, including the OS/PIE as shown here. De¬ 
spite the fact that the OS/PIE approach cannot reach (but will 
approach) the minimum of Q using noisy data because of 
the inherent properties of the algorithm, an OS/PIE-based re¬ 
construction strategy was adopted in the work presented in 
this paper because it provides band-limited reconstructions 

iniEi]. 
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FIG. 8. Suppl (S4). 3DBPP in the presence of noise, (a) Evolution of the cost function Q and of the error metric in the object space 
\\p — porigW"^ (b) shown as a function of iteration number for both the CG and the OS/PIE. (c) The cross-section through the center of the 
starting object porig and scattering geometry are shown, (d) Cross sections of the reconstruction p are shown at different iterations the CG and 
OS/PIE. 


























